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CALCTILATING  BODIES  T N SUPERSONIC  FLOW  AT  GREAT  ANGLES  OF  ATTACK 

Yu.  N.  D'yakonov,  L.  V.  Pchelkina,  I.  D.  Sandomirs kay a 

Tn  studying  rather  long  blunt  bodies  in  a flow  at  great  angles 
of  attack  a p (0  is  the  angle  of  slope  of  the  contour)  using 
various  numerical  Methods  adapted  for  calculating  smooth  gas  flows,  a 
feature  was  discovered  in  the  leeward  region  of  the  flow  f 1-3]  which 
made  it  impossible  without  taking  special  measures  to  continue 
calculating  for  the  quantity  x greater  than  a certain  critical  length 
vain*  *mp  (generally  depending  on  the  parameters  of  the  impinging  flow 
and  the  geometry  of  the  body).  Analysis  of  existing  results  and, 
foremost,  the  nature  of  the  distribution  of  the  peripheral  component 
of  velocity  vector  w close  to  the  surface  of  the  body  permits  us  to 
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assume  that  this  phenomenon  is  related  to  the  development  in  the 
shadow  region  of  the  flow  of  a feature  such  as  a shock  wave  resulting 
from  the  collision  of  supersonic  flows  directed  toward  one  another  at 
great  anqles. 

The  results  of  calculating  flows  near  sharpened  conical  bodies 
at  great  angles  of  attack  [4]  also  indicate  the  possibility  of  the 
development  of  such  a peculiarity. 

During  experimental  study  of  the  flow  past  sharp  cones  and 
rather  long  blunt  cones  at  great  angles  of  attack  under  actual 
conditions  in  the  presence  of  viscous  forces  on  the  shadow  side  of 
the  body,  developed  separation  zones  are  formed  f *3  and  elsewhere!.  It 
is  interesting  to  note  that  in  the  peripheral  direction  the  flow 
pattern  greatly  resembles  a supersonic  gas  flow  past  a round 
cylinder.  Specifically,  the  value  of  meridional  angle  at  which  the 
separation  zone  is  formed  in  the  three-dimensional  case  is  not  in 
good  agreement  with  corresponding  results  for  the  two-dimensional 
flow.  Thus,  in  the  first  approximation  we  can  assume  that  even  the 
experimental  results  do  not  contradict  the  assumption  of  the 
development  of  a shock  wave  type  of  peculiarity  when  a ^ 3 on  the 
leeward  side  of  rather  long  bodies. 

At  great,  angles  of  attack  the  pressure  on  the  attack  side  of  the 
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body  is  generally  significantly  greater  than  the  pressure  or  the 
shadow  side  of  the  flow.  Moreover,  in  the  presence  of  develope 
separation  zones,  calculation  for  an  ideal  fluid  in  this  region  can 
differ  markedly  from  the  real  values  of  the  gas  dynamic  parameters 
because  of  the  strong  influence  of  viscous  effects. 

The  situation  may  prove  analogous  to  that  which  occurs  in  a 
supersonic  flow  past  a sphere.  Using  inviscid  gas  equations  we  can 
calculate  this  flow  up  to  the  value  of  0 — 160°  (polar  angle  fl  is 
reckoned  from  tha  critical  point).  However,  when  O>0OtP  — 110-120° 
the  obtained  solution  does  not  reflect  the  real  flow  pattern,  since 
this  entire  region  lies  within  the  flow  separation  zone.  Therefore,  a 
ractical  solution  to  the  problem  of  calculating  bodies  in  a 
supersonic  gas  flow  at  great  angles  of  attack  aiqht  be  tha  following. 

Outside  of  the  separation  zone  the  flow  is  calculated  using  the 
model  of  an  ideal  gas,  which,  as  we  know,  in  this  case  describes  the 
process  very  well.  The  distribution  of  pressure  in  the  separation 
zone  is  assigned  on  the  basis  of  experimental  data  and  obtained 
calculation  results  for  the  smooth  flow  region.  Since  pressure  in 
this  zone  is  significantly  lower  than  the  pressure  on  the  attack  side 
of  the  flow  and  as  analysis  of  calculation  results  shows,  even  a 
rather  rough  distribution  of  gas  parameters  in  the  separation  region 
will  not  lead  to  great  error  in  the  values  of  the  sun  aerodynamic 
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characteristics  of  hlunt  cones.  Thus,  we  can  assume  that  such  a 
system  will  provile  an  accuracy  in  calculating  the  flow  past  bodies 
at  great  angles  of  attack  which  is  sufficient  from  the  practical 
standpoint . 

In  this  stuly  to  calculate  a three-d iaensional  steady  gas  flow 
in  the  supersonic  range  we  used  the  grid  method  in  the  form  discussed 
in  f6].  Let  ,,s  loo*  at  soie  modifications  of  the  method  which  must  be 
made  in  order  to  perform  a calculation  according  to  the  system 
described  above  (symbols  of  [ 6 ] used  here). 

Le*  us  assuma  that  we  have  a uniform  supersonic  flow  against  an 
ai isyaaet r ical  body  at  angle  of  attack  a and  (x,  j,  >|>)  a cylindrical 
coordinate  system  bound  to  the  body;  the  x axis  coincides  with  the 
axis  of  symmetry.  Meridional  angle  is  reckoned  from  the  plane  of 
the  angle  of  attack  «,  while  the  value  ^ = 0°  is  assigned  to  the 

attack  side.  In  view  of  the  symmetry  with  respect  to  variable  it  is 
sufficient  to  examine  the  gas  flow  in  the  range  180°. 

Let  us  introdace  into  the  region  where  we  seek  the  solution  x > 
x0 , 0 ^ F ( 1,  0 4 ^ ( »#  • grid  with  nodes  (rn,  mh,,  Z h?)  , where 
x = Ax;  At  =•  AJ;  A,  = A\|>  - const;  n = nt,  n,  + 1,  n,  | 2,  . . . ; 

m = 0,  l,  2 M;  f - 0,  I,  2 A,/M  - |;  lhL  - n. 


The  system  of  difference  equations  is  written  for  the  point 
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[(,+t),:  ("  + t)*‘;  '*•]• 

For  approximating  derivatives  with  respect  to  ^ we  use  difference 
relationships  analogous  to  the  formulas  of  (6)  fro*  [§]: 

j_ 

t 

i 

4"  p (XJ.  + I,  /+ 1 — Xi+I.  I-I  -f-  Xm,/+|  — — l)  l>  ^ ^ ^ 

where  X = [ u,  v,  w,  v,  £ j represents  the  vector  column  of  unknown 
functions,  u,  v,  w - projections  of  the  velocity  vector  onto  axis  x, 
v,  and  \j>  , respectively;  <*£.0^-0,  a*  0 = 1. 
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For  simplicity  let  us  assuse  that  the  separation  region 
coincides  with  on?  of  the  beass  ♦/“M*  and  l*^L . Let  us  assuae,  for 
exaaple,  that,  this  beam  l - L - 1,  and  then  it  will  be  the  last  even 
bean.  No  matter  how  we  write  derivatives  with  respect  to  'P,  as  long 
as  we  remain  within  the  framework  of  the  explicit  iteration  systea 
for  that  variable,  the  main  calculation  algorithm  remains  virtually 
unchanged,  although  calculation  of  the  right  part  of  the  difference 
equations  will  be  done  soaevhat  differently.  For  approxiaa tion  of  the 
derivatives  with  respect  to  ♦ on  this  last  bean  ue  night,  for 

example,  assume  the  following  relationships: 


DOC  = 1 PO  1 


PAGE  6 


f N * 5=  — ! — {a  (Xm+'i,  | — Xm+'t.  I— t + Xmti  Xit (-l)  i~ 

Jm+±',  tot 

-f-  P (Xm  + I,  I — Xm-H,  l-l  "t"  X,,,./  Xm.i_|)li 

(J*-\n+  ^ L_  ra (3XiVi.  i - 4A"mV,. + X;V..  »-*  + 

Wf  L+J-.i  «*,  . 

2 

+ 3x:t!-4x:t/-,-f  x;+!-s)  + 

+ p(3Xi+,.i-4X:+i.«_i  +Xi+i,<-|-f-  3Xi./  4X3,./-,  + 


Formulas  (2)  and  (1)  are  obtained  by  means  of  "unilat  eral" 
difference  relationships,  and  when  a = 0 and  the  number  of  iterations 
0^2,  they  assure  the  first  and  second  order  of  accuracy, 
respectively.  It  is  assumed  that  the  gas  flow  in  region  ♦ >♦<  has  no 
effect  on  the  purely  gas  dynamic  reqion  i.e.,  no  additional 

conditions  need  be  imposed  on  the  last  even  bean. 


Results  from  experiments  on  and  calculations  of  a supersonic 
flow  past  bodies  at  great  angles  of  attack  in  many  cases  confirm  this 
assumption.  Specifically,  in  [1,  1,  7 and  elsewhere]  it  was 
discovered  that  at  high  values  of  meridional  angle  >l>  on  the  shadow 
side  of  the  flow  the  peripheral  components  of  the  velocity  vector  w 
can  reach  supersonic  values. 


If  in  the  last  even  bean  peripheral  velocity  everywhere  is 
greater  than  the  local  speed  of  sound,  then  the  proposed  system  is 
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essentially  base!  on  the  same  ideas  as  the  method  developed  in  [ 4 "J 
for  calculating  two-dimensional  steady  gas  flows  near  conical  bodies 
at  great  argles  of  attack.  In  this  case  such  a formulation  of  the 
problem  is  obviously  correct.  As  indicated  by  the  exoerienced  gained 
from  calculating  and  comparing  results  obtained  by  the  proposed 
system  and  the  standard  method  [6")  (in  the  region  where  data  can 
still  be  obtained  using  the  standard  difference  system),  this  system 
can  be  successfully  used  even  when  the  condition  indicated  above  is 
destroyed.  The  explanation  for  this  is  as  follows:  although  the 
peripheral  velocity  component  on  beam  is  also  lower  than  the 

local  speed  of  sound,  for  the  studied  flows  it  is  comparable  to  it  in 
maqnitode.  Therefore,  little  information  is  transmitted  from  the 
region  > M1/  to  thm  region  , and  it  has  virtually  no  effect  on 

calculation  when 

The  fact  that  we  can  determine  the  numerical  solution  (which 
corresponds  to  reality)  in  the  region  of  the  gas  flow  bounded  by  a 
certain  line  on  which  the  velocity  component  which  is  normal  to  it  is 
below  the  speed  of  sound,  but  comparable  to  it  in  magnitude  (despite 
the  fact  that  the  formulation  of  the  problem  in  this  case  becomes 
incorrect,  strictly  speaking),  is  also  confirmed  by  the  special 
methodological  calculation  conducted  in  [8  1.  Studied  in  this  last 
work  was  the  flow  pattern  which  develops  when  a supersonic  freely 
expanding  stream  strikes  a surface.  The  characteristic  feature  of 
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such  a flow  is  th»  formation  of  an  open  narrow  subsonic  flow  zone 
adjoining  the  shock  wave  in  front  of  the  obstacle.  Let  us  note 
certain  peculiarities  in  the  flow  past  blunt  bodies  at  high  anqles  of 
attack.  If  we  coapare  the  data  in  Piqs.  1-4,  where  we  see  the  flow 
patterns  for  a blunt  cone  in  an  equilibrium  air  flow  at  nuaber  = 
?0,  pressure  p^  = 0.001  atm  phys.  and  temperature  = 250°K,  and  the 
pressure  distribution  on  the  surface  of  the  cone,  we  see  a marked 
shortening  in  the  region  effected  by  the  blunting  as  the  angle  of 
attack  qrows.  when  a = 30°  (Fig.  4)  at  great  distances  from  the 
leading  edge,  the  nature  of  the  pressure  distribution  on  the  body  Ps 
with  respect  to  the  meridional  angle  ^ is  in  good  agreement  with 
corresponding  results  from  [4]  for  sharp  cones  at  great  anqles  of 
attack. 

The  use  of  tie  system  for  the  shadow  side  of  the  flow  virtually 
eliminates  the  significant  limitations  on  the  length  of  a calculated 
body  flying  at  a great  angle  of  attack,  which  existed  previously  (for 
example,  the  standard  method  of  calculating  a flow  past  the  same 
cone,  even  for  a = 10°,  not  to  mention  qreater  angles,  could  only  be 
done  up  to  x — 10)  . 

Since  on  the  attack  side  the  flow,  even  at  short  distances  from 
the  blun*  portion,  becomes  close  to  conical  (with  the  exception  of 
the  thin  vortex  layer  near  the  body),  there  develops  a zone  of 
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drastic  rest ructur  inq  from  the  "almost,  clinical"  flow  to  an  entirely 
different  type  of  flow,  corresponding  to  the  flow  past.  *he  leading 
portion  of  a body.  In  the  vicinity  of  this  zone,  in  the  field  and  on 
the  shock  wave,  flow  peculiarities  of  a different  tyoe  may  develop, 
to  study  them  we  must  conduct  calculations  with  a dense  grid  for 
coordinates  5 and  along  with  a thorough  analysis  of  the  flow  in  the 
compressed  layer  and  the  geometry  of  the  shock  wave. 

Naturally  the  obtained  results  correspond  to  the  actual  flow 
pattern  only  in  the  gas  dynamic  region  inside  the  separation  zon». 
Based  on  the  data  of  experiments  in  [5]  for  cones  with  a half-angle 
of  3 — 10°  at  Nach  numbers  of  fl  ^ ■=  6 and  great  angles  of  attack  th» 
separation  zone  lies  in  the  region  of  meridional  angle  values  'I’ 

1 4 0°. 


The  proposed  calculation  system  can  be  easily  extended  to  the 
case  where  the  boundary  of  the  separation  zone  is  assigned  in  the 
form  of  a certain  function  (determined,  let  us  say,  from  the  results 
of  experimental  studies)  or  is  found  by  the  gasdynamic  parameters 
obtained  in  the  process  of  solving  the  problem.  This  is  especially 
true  if  we  consider  the  fact  that,  in  view  of  the  above  discussion, 
it  is  not  necessary  to  assign  this  boundary  with  a high  degree  of 
accuracy. 
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Further  refinement  of  the  formulation  of  the  problem  and 
obtained  results  requires  that  we  conduct  systematic  calru  lat.ional 
investigations  and  a thorough  experimental  study  of  the  shape  of  the 
boundaries  of  the  separation  zone,  distribution  of  qas  parameters  in 
♦■his  region,  and  the  effect  of  bodies  flying  at  great  angles  of 
attack,  the  parameters  of  the  impinging  flow.  Re  numbers,  the  nature 
of  the  gas  flow  in  the  boundary  layer,  etc.  on  the  flow  pattern  on 
the  shadow  side. 
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